A Biophysical Model of the Mitochondrial Respiratory System and Oxidative Phosphorylation

A computational model for the mitochondrial respiratory chain that appropriately balances mass, charge, and free energy transduction is introduced and analyzed based on a previously published set of data measured on isolated cardiac mitochondria. The basic components included in the model are the reactions at complexes I, III, and IV of the electron transport system, ATP synthesis at F(1)F(0) ATPase, substrate transporters including adenine nucleotide translocase and the phosphate–hydrogen co-transporter, and cation fluxes across the inner membrane including fluxes through the K(+)/H(+) antiporter and passive H(+) and K(+) permeation. Estimation of 16 adjustable parameter values is based on fitting model simulations to nine independent data curves. The identified model is further validated by comparison to additional datasets measured from mitochondria isolated from rat heart and liver and observed at low oxygen concentration. To obtain reasonable fits to the available data, it is necessary to incorporate inorganic-phosphate-dependent activation of the dehydrogenase activity and the electron transport system. Specifically, it is shown that a model incorporating phosphate-dependent activation of complex III is able to reasonably reproduce the observed data. The resulting validated and verified model provides a foundation for building larger and more complex systems models and investigating complex physiological and pathophysiological interactions in cardiac energetics

Modeling of Oxygen Transport and Cellular Energetics Explains Observations on In Vivo Cardiac Energy Metabolism

Observations on the relationship between cardiac work rate and the levels of energy metabolites adenosine triphosphate (ATP), adenosine diphosphate (ADP), and phosphocreatine (CrP) have not been satisfactorily explained by theoretical models of cardiac energy metabolism. Specifically, the in vivo stability of ATP, ADP, and CrP levels in response to changes in work and respiratory rate has eluded explanation. Here a previously developed model of mitochondrial oxidative phosphorylation, which was developed based on data obtained from isolated cardiac mitochondria, is integrated with a spatially distributed model of oxygen transport in the myocardium to analyze data obtained from several laboratories over the past two decades. The model includes the components of the respiratory chain, the F(0)F(1)-ATPase, adenine nucleotide translocase, and the mitochondrial phosphate transporter at the mitochondrial level; adenylate kinase, creatine kinase, and ATP consumption in the cytoplasm; and oxygen transport between capillaries, interstitial fluid, and cardiomyocytes. The integrated model is able to reproduce experimental observations on ATP, ADP, CrP, and inorganic phosphate levels in canine hearts over a range of workload and during coronary hypoperfusion and predicts that cytoplasmic inorganic phosphate level is a key regulator of the rate of mitochondrial respiration at workloads for which the rate of cardiac oxygen consumption is less than or equal to approximately 12 μmol per minute per gram of tissue. At work rates corresponding to oxygen consumption higher than 12 μmol min(−1) g(−1), model predictions deviate from the experimental data, indicating that at high work rates, additional regulatory mechanisms that are not currently incorporated into the model may be important. Nevertheless, the integrated model explains metabolite levels observed at low to moderate workloads and the changes in metabolite levels and tissue oxygenation observed during graded hypoperfusion. These findings suggest that the observed stability of energy metabolites emerges as a property of a properly constructed model of cardiac substrate transport and mitochondrial metabolism. In addition, the validated model provides quantitative predictions of changes in phosphate metabolites during cardiac ischemia

CrossTalk opposing view: Guyton's venous return curves should not be taught

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/101776/1/jphysiol.2013.260034.pd

Specification, Construction, and Exact Reduction of State Transition System Models of Biochemical Processes

Biochemical reaction systems may be viewed as discrete event processes characterized by a number of states and state transitions. These systems may be modeled as state transition systems with transitions representing individual reaction events. Since they often involve a large number of interactions, it can be difficult to construct such a model for a system, and since the resulting state-level model can involve a huge number of states, model analysis can be difficult or impossible. Here, we describe methods for the high-level specification of a system using hypergraphs, for the automated generation of a state-level model from a high-level model, and for the exact reduction of a state-level model using information from the high-level model. Exact reduction is achieved through the automated application of symmetry reduction and invariant manifold reduction techniques to the high-level model, allowing potentially significant reductions without the need to generate a full model. The application of the method to biochemical reaction systems is illustrated by models describing a hypothetical ion-channel at several levels of complexity. The method allows for the reduction of the otherwise intractable example models to a manageable size

A pH-Dependent Kinetic Model of Dihydrolipoamide Dehydrogenase from Multiple Organisms

AbstractDihydrolipoamide dehydrogenase is a flavoenzyme that reversibly catalyzes the oxidation of reduced lipoyl substrates with the reduction of NAD+ to NADH. In vivo, the dihydrolipoamide dehydrogenase component (E3) is associated with the pyruvate, α-ketoglutarate, and glycine dehydrogenase complexes. The pyruvate dehydrogenase (PDH) complex connects the glycolytic flux to the tricarboxylic acid cycle and is central to the regulation of primary metabolism. Regulation of PDH via regulation of the E3 component by the NAD+/NADH ratio represents one of the important physiological control mechanisms of PDH activity. Furthermore, previous experiments with the isolated E3 component have demonstrated the importance of pH in dictating NAD+/NADH ratio effects on enzymatic activity. Here, we show that a three-state mechanism that represents the major redox states of the enzyme and includes a detailed representation of the active-site chemistry constrained by both equilibrium and thermodynamic loop constraints can be used to model regulatory NAD+/NADH ratio and pH effects demonstrated in progress-curve and initial-velocity data sets from rat, human, Escherichia coli, and spinach enzymes. Global fitting of the model provides stable predictions to the steady-state distributions of enzyme redox states as a function of lipoamide/dihydrolipoamide, NAD+/NADH, and pH. These distributions were calculated using physiological NAD+/NADH ratios representative of the diverse organismal sources of E3 analyzed in this study. This mechanistically detailed, thermodynamically constrained, pH-dependent model of E3 provides a stable platform on which to accurately model multicomponent enzyme complexes that implement E3 from a variety of organisms